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TECHNICAL NOTE 2955 


ESTIMATION OF FORCES AND MOMENTS DUE TO ROLLING FOR 
SEVERAL SLENDER-TAIL CONFIGURATIONS 
AT SUPERSONIC SPEEDS 


By Percy J. Bobbitt and Frank S. Malvestuto, Jr. 
SUMMARY 


The velocity potentials, span loadings, and corresponding force 
and moment derivatives have been theoretically evaluated for a number 
of slender-tail arrangements performing a steady rolling motion at 
supersonic speeds. 


The method of analysis is based upon an application of conformal- 
transformation techniques. The utilization of these techniques allows 
the simple determination of the complex potentials for various types 
of two-dimensional boundary-value problems. 


In addition, two simple and often-used approximations to the 
rolling derivatives have been compared with the corresponding exact 
values determined by the method presented in this report. 


In order to show the importance of wing-tail interference, the 
effect of the flow field behind a rolling wing on the tail character- 
istics has been illustrated for a simple wing-tail arrangement. 


INTRODUCTION 


A definite need has been indicated for additional information on 
the contribution of various tail configurations to the lateral dynamic 
stability of airplanes and missiles at supersonic speeds. This need 
is especially acute with regard to information on tails which have only 
one plane of symmetry. Inasmuch as tail arrangements, in general, may 
be classed as nonplanar systems, the problem of estimating theoretically 
the aerodynamic loading for a prescribed motion, such as rolling, may 
entail some difficulty, particularly if the panels comprising the tail 
are broad and have subsonic leading edges. 
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The difficulty of obtaining rigorous solutions for three-dimensional 
nonplanar systems led to the slender-body theory. This approach con- 
siderably simplifies the analysis of such systems by allowing the neg- 
lect of compressibility effects and, as shown in a number of papers 
(refs. 1 to 5), permits a solution to be obtained by an evaluation of 
the incompressible disturbance potential in crossflow planes, which 
are defined as planes normal to the direction of relative flow. The 
application of slender-body theory has led to many worth-while results 
for a number of practical wing, wing-body, and tail-body configurations 
composed of pointed low-aspect-ratio lifting surfaces and slender 
pointed bodies (refs. 1 to 6). The damping-in-roll derivative for a 
slender-cruciform arrangement has been evaluated by Ribner in refer- 
ence Ó and Adams in reference 3. A velocity-potential solution reported 
by Westwater (ref. 7) is used in reference 6 to give the damping-in- 
roll solution for a slender configuration consisting of an arbltrary 
number of symmetrically placed triangular panels intersecting in a 
common chord. 


The present paper has two purposes. The first purpose is to 
determine the rolling and yawing moments due to rolling for a number 
of practical slender-tail arrangements. These arrangements consist of 
triangular-plan-form panels which intersect in a common chord but are 
not necessarily symmetrical with respect to this chord. The second 
purpose is to present the method used in solving the problem of the "4 
rolling two-dimensional boundary that was reported by Bickley in refer- 
ence 8. This method, in which the conformal-transformation technique is 
used, allows the solution of a wide range of two-dimensional problems 
and is especially applicable when the contours of the tail arrangements 
in the crossflow planes are not symmetrical. Also included herein are 
the necessary transformations and the integral form of the potential for 
a tail arrangement consisting of slender triangular fins attached sym- 
metrically to a slender circular cylinder. 


In order to assess properly the contribution of the tail to the 
rolling and yawing moments due to rolling for a complete configuration, 
the effect of the induced flow behind a rolling triangular wing on the 
contribution of the tail is approximately determined. 


For the reader whose primary interest is in the results rather 
than in the method, illustrative variations of the rolling stability 
derivatives for several series of tail arrangements, as well as sample 
span loadings and pressure distributions, are included in the figures. 
The formulas from which these loadings and derivatives have been 
obtained are presented in the section devoted to applications. 


The greater part of this investigation was carried out early in 
T992 
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SYMBOLS 


The positive directions of the forces, moments, and velocities 
are shown in figure 1. 


Xy 12 


Cartesian coordinates (original space coordinates) 
complex variable, z + iy 


complex variable, E + ir 


coordinates in complex (-plane 


spanwise variable in plane of wing 


Y 2 
b/2 


nondimensional spanwise variable in plane of wing, 


nondimensional coordinates, yy eu Z1 dr 


perturbation velocity components in X-, y-, and z-directions, 


respectively 
magnitude of vector in Z-plane 
air density 
flight velocity 


local pressure difference across surface 
free-stream dynamic pressure, = eve 


pressure coefficient 
free-stream Mach number 
perturbation velocity potential 
stream function 

complex potential, g + İy 


angle of sidesiip 


H,H1,B 
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rate of roll, radians/sec 


ratio of span of vertical panel to span of horizontal tail 
for inverted T-tail 


function of R, 


1/2 
1 + Vl Ur? / E E + cos y 


2g? l - cos 64 


ratio of span of vertical tail to span of horizontal tail 
for cruciform tail 


function of K, (KÊ + i) S erdt ES 
l - cos 01 
root chord 
spanwise variable along an arbitrary panel 
span of an arbitrary panel 
span of upper vertical panel 
span of lower vertical panel 
span of horizontal tail 
dimensions of general cruciform tail obtained by transformation 
arbitrary length 


arbitrary area 


pe 


aspect ratio of horizontal tail, | ———————————————— 
Horizontal-tail area 


(Vertical-tail span)? 
Vertical-tail area 


aspect ratio of vertical tail, 


_ (Twice span of largest panel)? 


Twice area of largest panel 
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ro body radius 

E spanwise circulation 

91 angle between root chord and leading edge of vertical panel 
99 angle between root chord and leading edge of horizontal panel 
94,92 angles used in transformation of general cruciform from 


Z-plane to €-plane 


T angle that panel makes with x,y plane 
N ' normal force 

x lateral force 

E rolling moment about x-axis 

N yawing moment about z-axis 

ON normal-force coefficient 

Cy rolling-moment coefficient, L'/q82 


(2) > 
C = — Á —— —— Á— 
¿CT ab(Horizontal-tail area) 


qe 
C1 — 
( don qo(Horizontal-tail area) 
r* 
G = RASS 
(2) um q(Twice area of one panel)(Twice span of one panel) 


z L 
q(Twice area of largest panel)(Twice span of largest panel) 


O 
Na 
d 


C7 
cx 


8V/. 0 
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Cw 
Cto) 7 op (Span of one panel) p—»0 


C O ee Ce DU UDIN E EISE 
( lp ve ap (Span of largest panel) p—0 


Oy lateral-force coefficient, Y/aS 
Y 
Cy — —M — ——— 
Jer a(Vertical-tail area) 


ye 
a(Vertical-panel area) 


Oe aaa 


2q(Area of one panel) 


ler = - 


2q(Area of largest panel) 
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" i 3(Cy) CT 
( p)oT aP (Vertical-tail span) 
V 2 pro 
E AC) 
(ran 7 
d= (Vertical-panel span) TER, 
a(c 
c - ( Y), 
P/ VT AÑ 
(Span of one panel) RS 
XC)... 
(ero) = " 
(Span of largest panel) p—>0 
Cn yawing-moment coefficient, N/qs1 
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(Yal li [NR EN, 

Subscripts: 

TE trailing edge me 

LE; leading edge 

CT cruciform tail 

TT inverted T-teil 

VT V-tall 

pal Y-tail 

ZV real and imaginary parts of a function, respectively 


GENERAL ANALYSIS 


Linearized three-dimensional supersonic flow is governed by the 
classical partial-differential equation 


(M^ - 1). ~ Pyy - zz = 0 (1) 


where @ is the disturbance velocity potential. If the disturbance 
body is slender in the axial direction, the potential flow is, toa 
satisfactory approximation (see refs. 4 and 9), controlled by the two- 
dimensional Laplace equation 


E + Pp, =0 (2) 


Briefly, the classification of a body as slender means that the 
production of perturbated flow in any cross plane (y,z plane) normal 


AW 
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to the free-stream direction (x-axis) is dependent only upon the dis- 
turbance contour in that plane and, hence, independent of the flow 
action in the preceding cross planes. Such a consideration implies 
the neglect of the gradient of flow in the longitudinal direction and 
necessitates the removal of the (MÊ - 1)¢,, term from equation (1). 
It seems reasonable to expect that the application of slender-body 
theory should give with satisfactory accuracy the aerodynamic loading 
for slender configurations for which the boundary conditions (such as 
rolling) do not depend explicitly on the longitudinal distance. 


Boundary Conditions 


Sketches of the tail arrangements considered in this paper are 
presented in figure l. Figure 2 presents a cross-sectional view of 
these arrangements and the prescribed boundary conditions for steady 
rolling motion. In the crossflow planes represented in figure 2, the 
two-dimensional potential flow satisfies the Laplace equation (eq. (2)), 
and the perturbated v and w velocities external to the disturbance 
contours must be continuous except possibly at the edges of the 
contour. As infinity is approached, the potential and the v and 
w velocities must become zero. Across the contour the normal flow 
velocity is continuous; whereas, the tangential flow velocity is 
discontinuous. 


Determination of Velocity-Potential Function 


The determination of the velocity potential in the crossfLlow 
planes for the tail arrangements shown in figure 2 can be simplified 
by a conformal transformation of the regions external to the disturbance 
contours to an upper half plane, the closed disturbance contours being 
transformed into the infinite real axis of the transformed plane. The 
velocity potential is then readily obtained in this plane, and by an 
inverse of the transformation, the expression for the velocity potential 
is established in the original plane. A brief discussion of the 
required conformal-transformation theory with particular emphasis on 
a rotating boundary (such as a rolling tail) is presented in the 
following section. 


Conformal transformations and auxiliary functions.- Any one of 
the crossflow planes for the tail configurations shown in figure 2 is 


called the original plane and designated symbolically as the Z-plane. 
The t(-plane (€ =£ + in) is designated as the complex plane in which 
the real axis is the transformed contour of the original plane and 

the upper half plane is the transformed region exterior to the contour 
in the original plane. 
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As is well-known, the transformation of an arbitrary, closed, recti- 
linear polygon from the Z-plane into the upper half of the f-plane can be 
expressed by the following form of the Schwarz-Christoffel transformation: 


dZ _ (Complex constant )TI(¢ - 0) 7/% (3) 


as (as ie 


where the position of the pole has been taken at L =i. The symbol y 
denotes the exterior angles of the polygon and c denotes the points 
in the €-plane which correspond to the vertices of the polygon in the 
Z-plane. The symbol II denotes the fact that each corner of the 


polygon requires a factor (¢ - o)7/x in the formula for dzZ/dt. The 
use of this expression requires that the numerator have no zeros and no 
branch points in the upper half of the f-plane and, consequently, none 
in the lower half plane. Integration of equation (3) yields 


|. g(t) 
7 go" (4 ) 


where g(t) is specified by the following relationship: 


(Complex constant )TI(¢ - c) 7/7 = g'(t)(t* +1) - etg(t) 


The determination of the g(t) function may be very difficult for 
some configurations. This transformation procedure is illustrated in 
reference 8. 


For convenience, the points defining the boundary contour in the 


Z-plane are denoted by 2+ iy; hence, the conformal relation between 
the original and transformed boundary becomes 


E g(t) (5) 


where t represents points along the real axis (the transformed boundary) 
in the €-plane. The preceding transformation may be expressed as 


Li 
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l 8,(8) + 1g (£) (6) 
t^ 3 
or 
ASS, _EylE) 
E E2 + 1 


Note that g,(t) and (E) are real-valued functions of E. 
„(E gy (E 


The following relationships are extremely helpful in the trans- 
formation of the boundary conditions from the contour in the Z-plane 
to the transformed contour (the real axis) in the €-plane, particularly 
when the contour in the original plane is experiencing a steady rolling 


motion. On the boundary in the Z-plane, ZZ = |Z|? = r^ where r is 
the magnitude of the vector-point function from the origin to the 


boundary. In the €-plane, the relation ZZ = r^ transforms into 


Lh) | where h(t) is defined as 
2 
(t£? + 1) 
2 2 
n(e) = [eG] * leis) (7) 
On the boundary Z = 3 ; hence, from this relation and equation (4) 
t- + 1 
zz = r2 = EBW) (8) 
(24 1)" 


(It should be noted that eq. (8) applies only on the boundary.) 
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and 
n(t) = & eC = [e (67]? + e? (9) 


which are the useful and necessary expressions that relate the 
original and transformed boundary conditions in the Z- and C-planes, 
respectively. Points on the boundary in the Z-plane coincident with 
the z or the y coordinate axis will take on values in the (-plane 


given by either ss LST or [Ey ()]^. 


As an illustration of the use of the preceding relationships, con- 
sider the following example: Given a boundary in the Z-plane which is 
rolling with steady angular velocity p about an axis through the 
origin of the boundary and normal to the plane of the boundary. The 
distribution of the prescribed normal velocity along the boundary is 
equal to pr, and this distribution must be equal to the tangential 
derivative of the stream function along the boundary. Hence, if the 


complex potential function is W - ( + iy, then = pre = y along the 


boundary. Now r2 = |z|? and, from the preceding discussion, the 
boundary condition for rolling in the Z-plane = plz]? transforms by 


| 1 h(g) 
use of equations (8) and (9) into = p —~———— along the transformed 
(12 + 1) 


boundary (which is the infinite real axis) in the (-plane. 


Integral expression for the potential in the C-plane.- As indi- 


cated in the previous section, the contour in the original plane is 
transformed into the real axis of the (-plane, and the region external 
to the contour is transformed into the upper half of the €-plane. In 
this plane the complex velocity potential W(t) = g(t) + 1y(£) is 
determined so that it satisfies the prescribed distribution of normal 
velocity along the boundary of the real axis. This normal velocity 
distribution, as shown previously, manifests itself through the stream 
function y(t) along the boundary. The remaining problem is to deter- 
mine a function of L that is finite and single-valued in the upper 
half of the t-plane with the imaginary part tending to y(t) as € 
approaches the boundary of the real axis from above. Bickley (ref. 8) 
presents the following integral function of 0 which satisfies these 
requirements: 


IS (10) 
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In the appendix this integral solution presented by Bickley is shown to 
be & special case of the general solution of the Riemann-Hilbert problem 
for the half plene which is reported by Muskhelishvili (ref. 10). For 


h 
the case of the rotating boundary, we ) becomes pats) as pre- 


o(g2 + 1) 


viously noted. The following integral expression is then obtained as 
the complex potential in the f-plane for the case of a rotating boundary 
in the original plane (see ref. 8): 


co 


wt) = E he)  _de — (11) 


Expressions for Pressures, Forces, and Moments 


Once the velocity potential has been determined subject to the 
required boundary conditions, the expressions for the lifting pressure 
end the corresponding forces and moments are easily obtained. The 
lifting-pressure coefficient is given by 


Ap. 24d. 2 aieo QR oT Mo (12) 


where IT, the spanwise circulation, is defined as du - Yo, the poten- 


tial jump across the surface. The subscripts simply denote values on 
opposite sides of a surface. Equation (12) is consistent with the 
small-perturbation theory only if the magnitudes of the perturbation 
velocities are equal across the lifting surface. When the magnitudes 
of the perturbation velocities are not equal across the lifting surface, 
equation (12) should contain differences in the squares of the v and 
w velocities. The squared terms contribute terms which are linear 
functions of p to the stability derivatives; therefore, these terms 
will vanish because the derivatives are to be evaluated as p— 0. 


The expression for the normal force acting on a panel of the tail 
configuration is 


So T BS 3 
N = oy f d Ol dg dx (13) 
o '"VL.E. 9X 
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or, since Ip - 0, 


N' = k 
= pV [pp ds (14) 
O 


Reduction to coefficient form by dividing by qS yields 


2 fs ) 
N= vg J, Prp ds (15 


The integral expression for the moment of the normal force is readily 
obtained on the basis of previously published analyses which have shown 
that, for rolling triangular arrangements with pressure coefficients of 
the form psSof(s/86), the longitudinal location of the center of pres- 


sure is at the Por point. With this stipulation, the yawing moment 
about the apex of the tail system may be expressed as 


at 
u 
EM 


So 
CroV sin T [ Top ds (16) 
O 


and division by gqS2 gives 


in T O 
B E [ Tom de 1 
n= 5 rg : TE (17) 


The rolling moment about the longitudinal axis is given by 


BO 
Bum ov | Dres ds (18) 
O 
and nondimensionally by 
Bo 
C1 = e Prps ds (19) 
P p81* O 
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For certain configurations it may be expedient, in order to 
facilitate integration operations, to determine the rolling and yawing 
moments by a consideration of the potential expression in the -plane 
instead of the original Z-plane. The necessary expressions for the 
moments are easily derived as follows: The rolling moment in the 
original plane can be expressed as 


ae o Orel dz + y dy) 


or 


B = Y Herz a( r2) 


Integration by parts gives 


L' = E iea - L ov dañs (208) 


In the t-plane equation (20a) may be expressed as 


s ii nle) ag 


when the contour in the original plane is assumed to be symmetrical, at 
least with respect to one of the coordinate axes. This condition of 
Symmetry, together with the antisymmetrical nature of the potential, 
causes the disappearance of the first term on the right-hand side of 
equation (20a) which expresses the rolling moment in the Z-plane. 
Similarly, the moment of the normal force (referred to the apex) may 
also be more easily obtained if the integration is performed in the 
C-plane as follows: 


5 5 ds g ds 
N = C» sin T ov | — de i — de (21) 
f Cr a Er R d T 


where a'b' and c'd' are those portions of the t-axis corresponding to opposite sides of a 


panel in the original plane. The side force, of course, 1s simply N divided by 7 cr. 


APPLICATIONS 


The results of the preceding analysis are applied in this section to the evaluation of the 
velocity potentials and associated force and moment derivatives for several tail arrangements. 
In addition, a brief discussion of the effect of the wing flow field on the tail characteristics 
is included. 


Forces and Moments on Rolling Tails 


General cruciform tail arrangement (symmetrical with respect to one plane).- The transfor- 
mation of the space outside a cruciform with vertical panels of different spans and horizontal 
panels of equal spans to the upper half of the C-plane is obtained from reference 11 as 


JT 
e| - cos 94)(1 - cos 95) + 2C (coc 91 cos 05 - 1) + (1 + cos 84)(1 + cos 5) | de 
4 cR 
re + 1 
(22) 


The orientation of the cruciform in the Z-plane and the correspondence of points in the Z- 

and t-planes may be noted from figure 5. The dimensions of the cruciform tail (that is, the 
panel spans H, H1, and B/2), as transformed by equation (22), are determined by the selection 
of cos 01 and cos 0o in the following equations (ref. 11): 


£ 
H = 2| - cos 03 ) a - cos H " (23a) 


OT 


GG63 NL VOVN 
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Hy. ela + COS 04)(1 + cos STS (23b) 


= cos 61 - cos 02 (230) 


To [t2 


In order to facilitate the selection of the cosine terms for known ratios 


of H/By and ze, equations (23) may be written in the following form: 


v HC 
i dE REIR 


cos 81 = 


e 
e E Lo 


(2) ) 


MU i de. m . |. a m 
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It is evident from the transformation that Z will be either real 
or imaginary on the boundary. The auxiliary function h(g) will 
therefore be given by 


h(&) = ai [ea - cos 84)(1 - cos 85) + 2t*(cos 04 cos 05 - 1) + 
(1 + cos 91)(1 + cos 02) (25) 


where the negative sign applies over the portions of the boundary from 


the points to © ana ÓN) to (2) (see fig. 3(c)) and the positive 
sign applies over the remaining portions of the boundary. Substituting 
these h(t) functions into equation (11) and performing the indicated 
integration results in the following expression for the complex potential 
in the C-plane: 


OC( EF, -= - 
G = D ML NS E AH (tan? Fz - Fy - 
ee (£8 + 17 
n Fz + F Fz - F 
tnt fp TR + M TI pe 
(^ + 1163 hL + 2-5ÀN tl 
FÉK__—_ a aaMaMMaaaaMMMMMMMlMŘħŰ e =Ma 
2 2 
(£€ + 1) (iz a dl NVR rau. e) 
(26) 
where 


F4 = (1 - cos 07)(1 - cos 02) 


Fo = cos 01 cos 92 - 1 


l - cos 0, cos Oa 
(1 - cos 04)(1 - cos 95) 


— 
— 


n, E cos 91 - COB 95 
(1 - cos 804)(1 - cos 95) 


If, in equation (26), € dis expressed in terms of Z (obtained from eq. (22)), there results 
the complex potential in the original plane (the Z-plane) for a cruciform with the dimensions H, 
Hi, and B/2 (see eqs. 23). 


The velocity potential in the t-plane (real part of eq. (26)), together with the h(t) func- 
tions (eq. (25)) and the relation in equation (20b), ylelds the rolling moment and, hence, the 
damping-in-roll coefficient as 


(Crp) op pee — eee kran E + Fy - ten"i [5 *) tnh ae 
L 


a(F4F) IT (L + Fz + &) (Qs Fz - Fy) 


[s (ant BTR eut Rel 12) 4] ef LE 


1 + F3 + Fy 


(Equation continued on next page) 


GG6e Ni VOVN 


6T 


ES + E, Pz E 
» 5 77h E 
D As IT 


m Fz + ts (1 + Fz - EF, )? 
tan-l fg A) - 5 - BEES? + hy 2). Ea Ai i 
z - Fy p Do po ES c SE. E 
vA Tal H pl + FE, +F, 


16(Fo - F4 ) + 12 


Es - Fy L 9 VF +R 


—————À + tre m0) + 38] - 
l + Fs ~ Fy, l + F3 - Fy (1 + Fz A pL )° 


9 = Gy + WP, + oF? - oF," 
pur de cl ae 


yFz - F | 
pp: i [x - aant ¡ETE - tan? [F5 - Fi) S 


E E OF Fs + AF + 16 
1aF + Fr 5 ^ = iot Pi + 16F5 + 18 - 
2 (1 + F3 + Fr) l + Fs + Fy 


(Equation continued on next page) 


GG62 MI VOVN 


Fx - F 16(Fo - F 12 
dFs-Tk | 16, NN AN E 
1 + Fa - Fy (1 + Fx - EY 1 + Fs - Fy 


1 + Fz + Py L + F3 - Fy 


yF qb JF =E 
ers + | LA EL M CER S 


-OF,2 + IRA + OF, - 2FQ^ - BFo - 9 an 
tan”? y Es er 29 in AU S: RU NEN MNA + Va TT A F 


E, MA a 
Je(Fo - Fy) t 15 8 (Fp KD + + 16(Fo - Fi) + 12 . l6R, - HEJÉ + 


(1+ F, + NE 1 + Fz + Fy 


Fx - F 2(Fo - F 16 
E A NN. MEA. eee ate 


1+ 45 - Fy |, Ps - ui" (1 + Ps - sL 


dI mi E, Y + em, | T 16(Fo = F4) + 12 Xx - 4 (tant yEs + Fy, = tan”? U P< se Fi, ) 
A A A Mes MEUM E e NENNEN Mh LE 
U 


l + Fz - Fy 
(27) 
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Td 


| B 
The variation of (C25) oA (eq. (27)) with the ratio E for & fixed value of H/H] of 2 


is shown in figure +. 


The side force resulting from the rolling motion may be determined utilizing the velocity 
potential and equation (21). In coefficient form this side force 1s given by 


PS ~ Fh co 
AY Gr) + G(t)at (28) 
5 0) + F 
(VF, + RL + 22 A) L ERE 


rala 7 


where 


a(t) = MO ] =o - (tan) RTE, - tani ¡EE 


eri Ur aa) 


EE ES as y 


L+Fz+Fy lc Fs -Fy 


h pu (fF Fi + BS =P EDI Gi - F0 + ( - Fo VIL 
(tS + s (= + Wye ae Vis + Fy - e ee + 1) sit + 2p. te - F} + 2F5 + 4 


CG6¢ NL VOVN 
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The variation of (Lo on /Av with the ratio == for a fixed value 


of H/H, of 2 has been determined by a numerical integration in the 


original plane and is presented in figure 5. With the lateral force 
known, the stability derivative Ong of the cruciform tall is easily 


obtained as pointed out in the section dealing with the pressures, 
forces, and moments. 


Of special interest are the limiting expressions of equation (22) 
which represent the transformations to the €-plane of: (a) the cross- 
flow plane of a pair of mutually bisecting, perpendicular surfaces and 
(b) the crossflow plane of an inverted T-tail. 


Cruciform composed of two mutually bisecting, perpendicular, 
triangular surfaces.- In treating the cruciform with equal-span verti- 


cal panels, the cross section may be considered as having the same 
orientation in the complex plane as in the real plane; that is, Hi is 


set equal to H and the horizontal tail has a span of 2H in both 
planes. In order that Hi be equal to H, cos 0o must be set equal 


to -cos 04 in equations (23) with the result that 


E = Hi = 2/L - cos?64 


and 


= 2 cos 01 


PO [Ed 


If the same substitution is made in equation (22) and the local span 
of the horizontal panel is denoted as b/2, the transformation (eq. 22) 
may be written in the more useful form 


2 
2 e y 2(1 + cos^01) al 
e l - cos^0j 
Zo Se C CE EDEN (29) 
C- + 1 


The values of L at points (2) and ÓD (fig. 3(c)) are obtained from 


1/2 /2 

ht 006: 8 l - cos 9 

equation (29) as Boucle and | ————— , respectively. 
L = COB 963 i 4 0605 07 
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1/2 


l + cos 94 
Let the parameter | ——————— 


be defined as e. Equation (29) 
l - cos 94 


may then be simply expressed in the form 


t$ yt : (e + s m 


t= + 1 


Z = 


If the vertical panels are now given a span of K 2, then 


e = ne +1-K (51) 


and 
e? + takki (32) 
e2 


The velocity potential for this case may be determined from 
equation (26) as 


_ (o/o lee - 1) (e2 + 1) le 1 
oe n ET E 


Y 2 150 
C -[e^ + =) +1 E 
ia oi ah log Cre can (33) 


e 


(2 410) (1 - et)(e + €) 


Substituting for € in equation (33) the equivalent in terms of Z 
(determined from eq. (30)) yields the potential in the original plane 
for a cruciform of arbitrary dimerisions. In the limiting cases where 
K=O or K= 1, this expression for the velocity potential reduces 
to that for the potentials given in references 2 and % for the rolling 
triangular wing and rolling equi-panel-span cruciform, respectively. 


From equation (33) the pressure distributions and span loadings 
for K values of 0, 0.5, 0.75, and 1 have been determined and are 
presented in figures 6 and 7. 
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The damping-in-roll coefficient for this group of tails, which is 
obtained from equation (27) by making the proper substitutions for the 
cosine terms in terms of e, is found to be 


2 h _ 6e? 
e e dado oes 


n 
(tan-le - ban") Z + Da + e2) . CONES 


Y 
= (54) 


het 
The variation of (297. pa with K is plotted in figure 8. 


Inverted T-tail.- By setting cos 09 = -1, the dimension Hj of 
equation (23b) vanishes and the cruciform in the complex plane (fig. 3(b)) 
degenerates into the inverted T-tail. If the value of L at point (2) 
for this case is designated as a, equation (22) reduces to 


yo ty = al 


| Ves (x1 


7, (35) 


The proportions of the inverted T-tail will then be governed by the 
selection of a in the following equations: 


B Lal ( 
— E ae 368) 
e L + a^ 
T (566) 
1 + ac 
g 2 
R = H = ye (360) 


to 
© 
PO 
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Equation (36c) may also be written as 
1/2 
1 + yl + ARE ` 
& = |- (36d) 
ORS 


to allow for the selection of a fora desired value of R. 


The velocity potential for an inverted T-tail of dimensions H 


and B/2 is found by a reduction of equation (26) and may be expressed 
as 


2 | 160p — ier mu uie 
(1 + a^ )x (€? ^H ijt b-a 


2 2 _ 2 E 
Dru Rd LER TM tanla) E NN (37) 


Sid: E) Fed 


A more general potential expression which will be valid for tails of 
dimensions h and b/2 may be evolved by multiplying the right-hand 


side of equation (37) by (h/H)? or (v/B)* and by altering the 
transformation formula so that on the vertical panel 


1/2 
sb P E (a2 + 1) 
E (38a) 
2 
eE- (HF 
and on the horizontal tail 
E 2 1/2 
d d 2 
"x +2 t 2\/l =- |= 
Mn tete Gre) >- (era) (385) 


2 
ho. 2 =) 
a2 a2 + S 
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The spanwise pressure distribution may be determined from equa- 
tion (12) together with the potential for the rolling inverted T-tail 
of dimensions h and b/2. The pressure distributions for various 
values of R are plotted in figure 9; spanwise loadings are presented 
for these same values of R in figure 1O. 

The following expression for (Yal is obtained from equa~ 


tion (28) when cos 82 = -1: 


1/2 


7 a(5a + 3) (a2? + 1) (rz - 4 tenia) 
Qam O MA (39) 


SL a7 + T g 


Figure 11 represents the variation of (Xp) an f with the parameter R 


The damping-in-roll coefficient nondimensionalized with respect to 
the horizontal-tail area is given by 


No AGL + a)” [ (ah + 202 + 9)(x - de ane 
( Lp pe 8 E Ñ 


ena 


ala? - 9)(4 tanha - x) . (at + iha? + 9)a2 


(140) 
2 CDS 


The variation of (tp lan [^ with R has been shown in figure 12 to 


depict the effect the addition of a vertical panel to a slender tri- 
angular tail has on the rolling characteristics of the tail. 


Rolling V-tails.- The expression necessary to transform a V-tail 
and the region exterior to it to the infinite straight-line plane can 
be obtained from reference 11 and is 

l-g 
2 
| k(-t?) 


r= +1 


7, (41) 
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where c defines the "dihedral angle + (see fig. 13) as 


GE 


(42) 


and the span of the tail panels as 


a/2 


OL tJ 


ea eI + 3 (43) 


l -og 


The positioning of the V-tail in the complex plane, as well as the | 
correspondence of points in the Z- and €-planes, may be seen in figure 15. 


In order to note the effect on the forces and moments of changing 
the dihedral angle, tails with dihedral angles of 22.57, 45°, and 67.59, 
1n addition to the two limiting cases of 0% and 90°, have been investi- 
gated. The auxiliary function h(t) and the velocity potential in the 
t-plane are listed in the following table for each of the five tails: 


v2 2 - 36 + 3) +1 
(2.064) (t9 4 1y 


Z w(t + O45) +3 
aan 


-t(t* + 3) 
u(t ei) 
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In order to show the various end-plate effects, trailing-edge 
potentials, as well as span loadings, have been plotted in figure 14. 
For comparative purposes the span loading for the equi-panel-span 
cruciform has also been included. The variation of (Crp) vni and 


(Yea with T is plotted in figure 15. This figure shows that, 


theoretically, V-tails with dihedral angles in the immediate vicinity 
of 45° yield more damping than a cruciform tail which has twice the 
damping-surface area. 


Rolling Y-tail.- The velocity potential for a rolling Y-tail with 
any desired ratio of the span of the arms to the span of the stem 
(vertical panel) may be readily determined by the method presented in 
this report. The Y-tail considered in the following analysis has an 
angle of 90° between the two arms and an angle of 1359 between the 
arms and the stem (fig. 16(a)). 


The expression required to transform the Y-tail and the region 


exterior to it to the infinite straight-line plane can be obtained 
from reference 10 as 


1/4 
Ee - cos 01) | (t - 22)3/#(a¢) 1/2 (UC) 


LZ +1 


E 


1/2 

Loe Cos. Url 
where e, as before, is defined as | —————— . The orientation 
l - cos 04 
of the Y-tail in the Z-plane and the correspondence of the points in 
the Z- and C-planes are shown in figure 16. The dimensions of the 
Y-tail are obtained by the selection of cos 0, in the following 


equations: 


1/4 
= = Piven (1 + cos 01) (45a) 
1/4 
Hz a lea. - cos «Y / (45b ) 


20 NACA TN 2955 


The required h(t) functions are 


4 1/2(1 - cos Ele? - an 


and 


4 \/2(1 - cos o4) t (g2 - aj 


on the arms and stem, respectively. Substituting these functions into 
equation (20b) and integrating over the proper portions of the boundary 
(see fig. 16) yield the complex potential in the f-plane as 


aa Ny eee a A 
W=2 2(1 - COS 01)?» Cyl + e* (557 + 1 - 205) + 


2(t2 + 1) 
: (12 4 1)" (t2 e 1)? 


MIGR ae 


(ean): 


(46) 


The nondimensional.velocity-potential parameter on the arms is 


eta - COS 91) V1 + el 


d. 


p(b/2)* 5/5(1 + cos a )^ 


EE + 1 - 2e*) 


+ e 
2(t2 4 1)° Vi + e? 


-t* (4, L eH) 4 ef - 2 " tle? - (235/2 


e(t? + 1) V1 + ee(t2 + ijs 


(47) 
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and on the stem or vertical panel is 


(48) 


In the limiting cases cos 94 = l and -1 the potentials are obtained 


for the rolling V-tail with a dihedral angle of 15% and for a triangular 
panel rolling about one edge, respectively. Sample trailing-edge 
potentials and span loadings are presented in figures 17 and 18. 


For convenience, only numerical integrations were made of the span 
loadings to obtain the side force and rolling moment due to rolling. 
These results have been nondimensionalized and are presented as the 


parameters (C A' and (C A' in figure 19. 
to)yr Lp 


Figure 20 has been prepared to illustrate the relative damping 
qualities of the various tail arrangements of panels of the same span 
and area considered in this paper. The relative merit of the Y-tail 
as shown in this figure may be of particular interest. 


Cruciform tail-body arrangement.- The method presented herein can 
easily be applied to a cruciform tail-body arrangement to determine the 
integral solution of the complex potential in the C-plane. No effort 
was made to evaluate the integrals; therefore, the analytical expression 
for the potential and, hence, the associated pressures and span loadings 
are not presented for this case. Nevertheless, the development of the 
integral form of the complex potential in the €-plane was considered of 
sufficient interest to warrant presentation herein. 


Ihe contour of the cruciform tail-body configuration shown in 
figure 21(a) may be transformed to a pair of mutually bisecting, per- 
pendicular, triangular surfaces by an application of the Joukowski 
transformation (fig. 21(b)) and then to the infinite Straight line 
(fig. 21(c)) by use of a transformation similar to the one used in 


the previous illustrations. This procedure results in the following relation between Z and L: 


r 1/2 


(0/21 «rfl, 2 |k2(m/232 - rof] 
2/2 Ke [(o/2)2 + re] 
c2 + 1 


+1 ¿2 +1 


1/2 


(9/24 + roel” e zd 2 [Ke (9/27 ? zu + 
K^ Iv /2)° + r] 


te + L 


PEES lim ^ (p/2)*(t? + ij* 


On the vertical and horizontal panels both terms will be real or both imaginary; whereas, on 
the body the first term will be real and the second imaginary. 


As noted previously, in order to determine the complex potential, the h(t) function must 
be substituted into equation (11) and integrated over the boundary. The h(t) functions over 
the portions of the boundary corresponding to the horizontal and vertical panels (see fig. 21) 


will be plus and minus [5 (6) + to(£)]^, respectively, where f4(t) is the numerator of the 
first term of equation (49) and fo(t) is the numerator of the second term. The h(t) func- 
tion over the portions of the boundary corresponding to the body is [f COT? - AGA The 
integral equation for the complex potential is then 


(49) 


NON 
no 


GG6c NL VOVN 
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(gg + 16 - €) E (2 + 1) E - £) 


“e [e (8) + eale de IT 1E GJI R Ole, de I 


JO Ew «situ, [O (Bal - Ged fas, 


O (?+U%eE-t) “O (62 + 1) - £) 


ii [8(6) + fo(e)} at M a Ir. GT x Sat JI? dE 
Gaan tn “© (62 + 1) = t) 


JO Et «solu, [O [pa - feet?) as 


O (2+ e@-t) "O (62 + 158 - £) 


(50) 


lS [YLST + a 


@ (62+ 1% - £) 


Where the values of E 
(fig. 21(c)) are given by 


2 2 
E . € Ik? (0/2) - re? eH E " zal 


Ke (0/24 + r2]? 


2 
4 2 2 
B epeQeP - PD - 1 + 1600 + 
pie K^ [(b/2)* + rêl" 


at the points @) and @ by 


[Gra - «el 
E = £ |< = + 1 
re [(9/2) + ro2|* 


as determined from equation (49) at the points 


SMS" - ro? 


K R ye + ro? | 


and at the points (9) and (5) by the reciprocal of equation (52). 


BD, ©, ©, aa © 


1/2 
[x (0/2) : Lo] 
(5/2)? + 238 
(51) 
(52) 


1/2 
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Approximation of the Rolling Stability Derivatives 


For many arrangements, because of the inherent difficulties 
involved in making & rigorous estimate of the stability derivative CY, 


(ana hence, Clg and Enp} an approximation to this derivative based 


on a knowledge of the aerodynamic loading due to sideslip becomes 
necessary. In the simplest of these approximations, B is assumed to 
be equal to the average angle of attack due to rolling ph/2V. Ina 
second and more accurate approximation B is assumed to be equal to 
the spanwise distance to the center of pressure of the sideslipping 
tail multiplied by p/V. If Cy, is expressed in the form 


Cy, == = (23) 
P ^ og 263 
V 
the first approximation is 
1 
C S 0 
eee i (54) 
and the second approximation is 
Z 
C = =C 
Y 7705 Ya (55) 


where Z is the spanwise distance to the center of pressure. 


An indication of the accuracy of these approximations may be 
obtained by & consideration of the inverted T-tail. A theoretical 
estimate of Ye for an inverted T-tail may be determined from refer- 


ence 12 and can be expressed as 


2 
1 + 2a (56) 


X 
(roer y p^ l + ac 


This expression together with equation (39) provides an interesting 
theoretical check on the accuracy of the preceding approximations when 
applied to a slender inverted T-tail. 


36 NACA TN 2955 


The correction factor given by the ratio of equation (39) to 
equation (56) is 


(or - ha(589 + 5) + 3(82 + Dar - 4 tanla) 
1/2 


(T) 


(Ya Tn 6n(a? + 1)" (1 + 282) 


This ratio is compared in figure 22 with the approximations given by 
equations (54) and (55). The values of (CY¥p Tr obtained from equa- 


tions (54) and (55) are presented in figure 11. It may be noted that 
equation (55) gives excellent agreement (within 5 percent) with the 
exact values of (Yo) for values of R greater than 2. For values 


of R less than 2 an estimation of (Co) based on the approxima- 
tion obtained by equation (55) mey not be advisable. 


As a point of interest, the varlation of (CY alp /Av end (C26 Jony [^ 


(determined from ref. 12) with R is presented in figures 23 and 2}. 
The expression for (C1g)mr is 


C1 > d. Oe M 2(1 + a2 )* 2all - ac) meat 1 - al " «(l + a2)* 
TT = A A A ee EA 
( P i 5a alb (1 + a2)? l4 a? aÔ 


(58) 


Effect of the Wing Flow Field on the Tail Characteristics 


In order to make an accurate estimate of the derivatives ls 
and ois for a complete wing-tail-body configuration similar to 


present-day missiles, the effect of the wing flow field on the tail 
must be considered. The importance of these effects for a rolling 
wing-tail combination is simply illustrated by an approximate deter- 
mination of the downwash and sidewash velocities induced on a slender 
inverted T-tail by a slender, rolling, triangular wing (see fig. 25). 


Ihe perturbated downwash and sidewash velocities were obtained in 
& plane perpendicular to the free stream at an infinite distance down- 
stream of the wing. Although the wing is of low aspect ratio, the 
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trailing vortex sheet has been assumed to remain flat. Even though 

the vortex sheet may actually roll up, no theoretical or experimental 
evidence is available to indicate the aspect ratios or rotational speeds 
at which this rolling up will occur. A further assumption is that, for 
the usual values of pb/2V encountered and for tail lengths of prac- 
tical magnitude, the rotation of the vortex sheet at the tail location 
can be neglected. 


Based upon the preceding assumptions, the downwash and sidewash 
determined from the Biot-Savart equation are formulated as follows: 


HA 
dy 
v(o,y1,0) = Pg i ar- "EIC (59) 


]: dy 
Zi 1 ar 2 
v(9,0,21) = = aoe — —s— (60) 
b/2 -l dyo Z4 + yo 


For a slender triangular wing (ref. 5), the rate of change of the 
circulation with respect to yo is 


dyp TE 


Substituting equation (61) into equations (59) and (60), performing the 
integration, and taking the partial derivative with respect to pb/2V 
give 


ar L - 2y,£ 
—= 120 E (61) 


no 


= yi | Cyil «1) (62) 
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and 


Po 


p o 
a H-A - 2 (zi > 0) (63) 
ay l + ze 


The induced downwash (eq. (62)) will equal both in magnitude and direc- 
tion the downwash due to rolling of the horizontal tail which is 


pb y 
W = - = ese = V v c — 6l 
py Mi: (64) 


— E aH (65) 


The damping force resulting from the difference of these two perturbated 
velocities (eq. (65) minus eq. (62)) will therefore be zero. The induced 
sidewash from the wing exerts a side force on the vertical panel of the 
tall in opposition to the side force produced by the isolated tail; 
hence, the effect of the wing sidewash is to decrease the amount of 

tail damping. This decrease will be rather large for most wing-tail 
arrangements and it is possible for the induced force to be greater 

than the damping force of the isolated tail with the result that the 
vertical tail contributes an antidamping force to the system. For the 
illustrative example the induced side force on the vertical panel is 
approximately 82 percent as large as the side force due to rolling. 

The damping in roll of the tail immersed in the flow field behind the 
wing is of the order of 15 percent of the value predicted for the tail 
alone. The induced angle of attack on the vertical panel, the angle 

of attack due to rolling, and the resultant effective angle of attack 
are plotted in figure 25. 


CONCLUDING REMARKS 


A method, based on conformal-transformation techniques, for solving 
two-dimensional boundary-value problems has been used to evaluate the 
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velocity potentials, span loadings, pressure distributions, and asso- 
ciated stability derivatives for several slender-tail arrangements. 
Illustrative variations of the rolling stability derivatives for 
several series of tail shapes, as well as sample span loadings and 
pressure distributions, are included. 


Iwo simple approximations of the rolling stability derivatives 
have been examined in view of the exact values determined by the method 
presented in this report. In addition, the importance of wing-tail 


interference has been shown with the aid of some elementary flow-field 
calculations. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 17, 1953. 
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APPENDIX 
DERIVATION OF THE EQUATION FOR THE COMPLEX POTENTIAL 


The purpose of this analysis is to derive the equation (eq. (10) in 
the body of the paper) 


W(t) sgri- Ij BH 


from a knowledge of the general solution of the Riemann-Hilbert problem 
for the half plane. The Riemann-Hilbert problem is formulated as 
follows: A function W(t) = Q(&) + iy(6£) must be found which is 
regular in the upper half plane Ln » O), bounded at infinity, and 
satisfies, on the infinite real axis L, the boundary condition 


1(£)8(&) - mle)y(e) = ne) (Al) 


where 2(&), m(t), and n(t) are real continuous functions of £ 
given on L. The functions are assumed to satisfy the Holder condition 
on L (see refs. 10 and 13) and the additional condition that 

12 + m= 40 for all values of E including £ = o. In regard to the 
Holder condition, a function, say m(E), will satisfy this condition 

on L if, for any two points, Ey, Ep on L, 


lalea) - xe | S4] - 51] 


where A and p are positive constants. The constant A is called 
the Hólder constant and y the Holder index. This definition of the 
Holder condition has been obtained from reference 10. 


The general solution of the Riemann-Hilbert problem under the 
conditions stipulated has been presented by N. I. Muskhelishvili in 
the following form (see ref. 10): 


W(t) = A A (content x(t) (a2) 
^ x*() 16) + am(e)] (e - 6) 
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where X(t) is the fundamental solution of the associate homogeneous 
problem and 


SIR) = Hm X(0 


7—>0 


The fundamental solution X(t) in terms of the boundary conditions is 
given by 


where 
aor l- im 
L pues l + - 
(t) = i 
-90 E - Eo 
and 


EO = lim C 
M0 


Now, for the boundary conditions considered in this paper, it is 
necessary to set 1(t) =O and m(t) = -1; hence, n(e) = y(8). It 
follows that 


X(¢) = -1 


X(t) = -i 


and the general solution expressed by equation (A2) resolves into 


W(t) = = MICI + i(Constant) (A5) 


-ot- 


which is the Bickley relation (except for constent term) that is presented 
as equation (10) in the body of the paper. The constant term of equa- 
tion (A5), however, is zero for the boundary conditions considered herein, 


since the imaginary part of W(t) as LR must equal YLE). 
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Fredrick 


Figure 1.~ Typical tail arrangements considered in the analysis. Center 
tail arrangement illustrates the System of axes used and the positive 
direction of the velocities y Moments, and lateral force. 
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Figure 2.- Boundary conditions on tail arrangements for Steady rolling 
motion. 
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Real plane 


(a) Cross section of cruciform in physical plane. 
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Complex Z-plane 


(b) Crueiform oriented for two-dimensional analysis. 
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(c) Points in the C-plane corresponding to the vertices in the Z-plane. 


Figure 3.- Physical and transformed planes used in analysis of rolling 
tails. 
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Figure 4.- The damping in roll of a slender cruciform teil for various 
values of n and a fixed velue of 2 for the ratio H/H}. As 
1 


T — o, Cy a approaches the slender-tail value of -1/32. 
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Figure 5.- The effect on the side force due to rolling of a slender cruci- 
form tail of increasing the horizontal-tail span while maintaining the 
ratio of the upper-vertical-panel span to the lower-vertical-panel span 


equal to a. 
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(a) Horizontal panels. 


Figure 6.- Pressure distributions on the horizontal and vertical panels 
of four cruciform configurations with ratios of vertical-tail span 


to horizontal-tail span of 0, 0.5, 0.75, and l. 
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(b) Vertical panels. 


Figure 6.- Concluded. 
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(&) Horizontal panels. 


Figure /.- The span loading in roll on the vertical and horizontal panels 
of four cruciform configurations with values of K of O, 0.5, 0.75, 
and l. Negative loading indicates damping forces directed in negative 
y-direction on upper vertical panel and in the positive y-direction 


on lower vertical panel. 
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(b) Vertical panels. 


Figure 7.- Concluded. 
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Figure 8.- The effect on the damping in roll of a slender cruciform tail 
of varying the ratio of the spans of the two mutually bisecting perpen- 
dicular wings. 
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Figure 9.- Pressure distributions on the vertical and horizontal panels 
of a series of inverted T-tails. 
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(b) Vertical panels. 


Figure 9.- Concluded. 
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(a) Horizontal panels. 


Figure 10.- The span loadings in roll on the vertical panel and horizontal 
tail of & series of inverted T-tails. 
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Figure 10.- Concluded. 
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Equation (54) 
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Equation (55) 


Equation (39) 
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Figure ll.- The variation of the stabillty-derivative parameter S 
» 


with R for a rolling inverted T-tail. 
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Figure 12.- Effect on damping in roll of adding & vertical panel to a 
slender triangular tail. 
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(a) Cross section of the V-tail in the physical plane. 


Complex Z-plane 


(b) V-tail oriented for two-dimensional analysis. 
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(c) Points in the t-plane corresponding to the vertices in 
the Z-plane. 
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Figure 13.- Physical and transformed planes used in analysis of 
rolling V-tail. 
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Figure 14.- The span load distributions and trailing-edge potentials on 
the rightthand panel of a series of V-talls in steady rolling. The 
span loading on a panel of an equi-panel-span cruciform is included 
for comparative purposes. 
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Fi e 15.- Variati f the stability derivatives and E 
gur etion o y (Cx) y ( bp) VI 


with the dihedral angle T for a series of V-tails. The circled 
points indicate theoretical values. 
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(a) Cross section of Y-tail in physical plane. 


Complex Z-plane 


(b) Y-tail oriented for two-dimensional analysis. 
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(c) Points in the -plane corresponding to the vertices in 
the Z-plane. 


Figure 16.- Physical and transformed planes used in analysis of 
rolling Y-tail. 
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Figure 17.- The span load distributions and trailing-edge potentials on 
the inclined right-hand panels of a series of Y-tails in steady 
rolling. It should be noted that the small sketches on each graph 


are intended only to depict the ratio of the panel spans HE 


B/2. 
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Figure 18 
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.- The span loading in roll on the vertical panel of a series 
of Y-tails. 
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Figure 19.- The effect on the stability-derivative parameters (Or. loa 


and (DAT of changing the span of the vertical panel on a slender 


Y-tail. The derivatives are nondimensionalized with respect to the 
diameter of the swept-out circle and twice the area of the largest 


panel. 
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Figure 20.- Variation of s with aspect ratio for five slender config- 


urations. All panels of the configurations are of equal span and area. 
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(a) Crossflow plane of cruciform-tail——body combination 
oriented for twoedimensional analysis. 


Intermediate plane 


(b) Cruciform-tail—body combination transformed by Joukowski 
transformation. 
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(c) Points in the C-plane corresponding to the vertices in the Z-plane. 


Figure 21.- Physical and transformed planes used in analysis of rolling 
tail-body combination. 


NACA TN 2955 69 


All three curves approach 


. Equation (54) e this value at R = 9 
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Figure 22.- The variation of the ratio CENT with R for 


the inverted T-tail as determined by exact and approximate methods. 
At R = «wo all equations give ratio equal to 0.5. 


Figure 25.- The variation of Gar. with 
R for a sideslipping inverted T-tall. 
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Figure 24.- The variation of GA. with 
R for a sideslipping inverted T-tail. 
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Figure 25.- The effect of sidewash from a slender rolling wing on the 
angle-of-atteck distribution along the vertical panel of a T-tail. 
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